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Abstract 

Transient homogeneous nucleation is studied in the hmit of large critical sizes. Starting from pure 
monomers, three eras of transient nucleation are characterized in the classic Becker-D5ring kinetic 
equations with two different models of discrete diffusivity: the classic Turnbull-Fisher formula and 
an expression describing thermally driven growth of the nucleus. The latter diffusivity yields time 
lags for nucleation which are much closer to values measured in experiments with disilicate glasses. 
After an initial stage in which the number of monomers decreases, many clusters of small size are 
produced and a continuous size distribution is created. During the second era, nucleii are increasing 
steadily in size in such a way that their distribution appears as a wave front advancing towards the 
critical size for steady nucleation. The nucleation rate at critical size is negligible during this era. 
After the wave front reaches critical size, it ignites the creation of supercritical clusters at a rate 
that increases monotonically until its steady value is reached. Analytical formulas for the transient 
nucleation rate and the time lag are obtained that improve classical ones and compare very well 
with direct numerical solutions. 

PACS numbers: 82.70.Uv, 83.80.Qr, 05.40.-a, 05.20.Dd 
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I. INTRODUCTION 

Homogeneous nucleation occurs in many examples of first order phase transitions |]|] 
such as condensation of hquid droplets from a supersaturated vapor, glass-to-crystal trans- 



formations [2], crystal nucleation in undercooled liquids [3|, and in polymers colloidal 
crystallization j^, growth of spherical aggregates beyond the critical micelle concentration 
(CMC) jsl Tj^and the segregation by coarsening of binary alloys quenched into the misci- 




bility gap j8|, |9|, llO|- In condensed systems, a long time elapses before the nucleation rate 
(at which stable nucleii larger than the critical size are generated) reaches a steady state, 
therefore these systems offer excellent opportunities to study time-dependent nucleation . 

Understanding the kinetics of nucleation and growth beyond the determination of the 
steady-state nucleation rate is a task of great importance and not yet completely accom- 
plished. For example, it is desirable to obtain a simple asymptotic description of the tran- 
sient until the steady-state nucleation stage sets in. Moreover, there is no clear distinction 
between nucleation and growth, and an unified theory of both processes does not exist 
despite recent attempts at bridging the gap between nucleation and late-stage coarsening 



theories 
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In this paper, we consider the problem of describing the approach to steady-state nu- 
cleation within the classical nucleation theory [3J. Thus our starting point is the Becker- 
Doring (BD) discrete kinetic model of nucleation and indefinite growth of a stable phase 
from a metastable state 0, Q, . The BD model contains two kinetic rate constants that 
are related to each other by assuming detailed balance. To complete the description of the 
BD equations (BDE), a model for one of the rate constants, usually a discrete diffusivity 
describing the rate at which a cluster loses one monomer, is needed. In the classical theory, 
the discrete diffusivity is given by the TurnbuU-Fisher (TF) expression which assumes that 
a monomer has to overcome an activation energy barrier for its transfer across the interface 
of a cluster. The TF discrete diffusivity is therefore proportional to the surface area of the 
cluster [3| • Other models are selected so as to yield the known expression for the adiabatic 
growth of a nucleus of critical size by either diffusion or by heat transfer. The discrete 
diffusivity of these later models is proportional to the cluster radius. 

No matter which discrete diffusivity is used, starting from an initial condition of pure 
monomers surpassing the CMC, we expect that cluster size increases and stable supercritical 
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nucleii are formed at a nucleation rate that will eventually become stationary at an expo- 
nentially small value. After the stationary nucleation has set in, the supercritical clusters 
continue growing, and the discrete diffusivity of the BDE can be ignored in the description 
of their growth, which is a pure convection in the space of cluster size. For precipitation 
processes, this will eventually result in late stage coarsening which we will not study in the 
present paper. 

The small parameter that informs our asymptotic analysis is e = k~^^^, where is the 
number of monomers in the critical nucleus, a 'large' quantity that ranges between 20 and 
1000 for common materials j^. Using other small parameters, such as the supersaturation, 
yields particular cases of our results (cf. Chapter 2 by Neu and Bonilla in Ref. Q]). The 
analysis of the BDE in the limit as e — > distinguishes three well defined stages or eras in the 
approach to the stationary nucleation rate. Starting from the initial state of pure monomers, 
a continuous distribution of cluster sizes is established at the monomer's expense during the 
first era. During the second era, the clusters grow to the critical size in such a way that 
their size distribution is a traveling wave front in size space. As this wave reaches the critical 
size, the formation of supercritical nucleii starts, nucleation is ignited, and the nucleation 
rate increases from zero to its stationary value during the third era. We have obtained two 
different expressions for the nucleation rate (which is of paramount importance to compare 
with experiments): (I) a general expression in terms of the instantaneous location of the 
wave front and its instantaneous width, which solve two given differential equations, and (II) 
a more explicit description of the nucleation rate in terms of the solution of the linearized 
wave front position with an origin of time at the time Im needed for the exact wave front 
to advance from pure monomers to a certain near critical size. Numerical solution of the 
model confirms all the theoretical predictions. 

Most previous studies of transient nucleation considered the Zeldovich-Frenkel equation 
(ZFE), which is a Fokker-Planck-type equation resulting from taking the continuum limit 
of the BDE[1X]- Zeldovich 17] set the discrete diffusivity equal to its value at the critical 
cluster size and used a parabolic approximation for the variation of the freejsnergy. The 
resulting expression for the transient nucleation rate was rather inaccurate Until the 
mid 1980s, work on the ZFE was based on similarly uncontrolled approximations Some 
of them gave expressions for the nucleation rate and time lag close to the values obtained 
by numerically solving the BDE for particular parameter values, but were far off for other 
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parameter ranges Asymptotic theories for the ZFE were elaborated later 
There are two main differences between asymptotic results obtained for the discrete BDE 
and those obtained for the continuum ZFE: (i) the time lags for transient nucleation are 
different, as explained by Wu [ij], and (ii) the width of the wave front and the time to 
ignition are different (wider for the ZFE). Nevertheless, other magnitudes such as relaxation 
times and the stationary nucleation rate are the same for asymptotic approximations of both. 
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yields expressions for the nucleation 



21[ and by Shi et al 



22j |. although their 
. For large critical 



the BDE and the ZFE. Thus our simplified theory 
rate that are similar to those found by Shneidman 
time lags differ from ours, as one would expect from Wu's arguments |1 
sizes, our approximation (I) is better. 

The rest of the paper is as follows. In Section IHl we review the Becker-Doring model 
for nucleation and growth of spherical aggregates with the Turnbull-Fisher (TF) discrete 
diffusivity The binding energy of the aggregate with k monomers {k cluster) relative 
to isolated monomers in solution is {k — 1) times the monomer-monomer bond energy plus 
a term proportional to the surface area of the aggregate. Beyond a critical density no 
equilibrium size distribution exists and the aggregates grow indefinitely. The main results of 
our asymptotic analysis are derived in Section UTTl and compared with the numerical solution 
of the BDE with the TF discrete diffusion coefficient describing devitrification of lithium 
disilicate glass. Our results compare favorably with previous theories based on the ZFE, the 
continuum approximation of the BDE. However, when compared with experimental data for 
glass disilicate, the theoretical time lag is about 30 times smaller. To improve the agreement 
with experiments, we propose in Section HVl a different discrete diffusion coefficient selected 
so as to yield the known expression for the adiabatic growth of a nucleus of critical size 
by heat transfer. The asymptotic theory for the resulting BDE is similar to that explained 
in Section IIIIl and the resulting time lag is much closer to experimental data. Section 
compares our asymptotic results for the transient nucleation rate and for the number of 
supercritical clusters to previously known analytical formulas (unfortunatelv all of them 



dealing with the continuum ZFE, not with the discrete BDE as ours do) [2d,|2l|,|2d,|2i|2|. 
Technical matters are relegated to Appendices. 
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II. KINETIC EQUATIONS AND STATIONARY SOLUTIONS 



The model presented here is nucleation in a lattice in which there are many more binding 
sites, M, than particles, A^, 0. We shall consider the thermodynamic limit, N ^ oo with 
fixed particle density per site, p = N/M. Let pk be the number of clusters with k particles 
or, in short, k clusters, and let pk = Pk/M he the density of k clusters. Note that the number 
densities per site, p and pk, are both dimensionless. Number densities per unit volume are 
obtained dividing p and pk by the molecular volume, v = V/M. Particle conservation implies 
that the total particle density p is constant: 

oo 

= p- (1) 

k=l 

In the Becker-Doring kinetic theory of nucleation, a k cluster can grow or decay by capturing 
or shedding one monomer at a time. Then j?! 

Pk = jk-i - jk = -D- jk, k> 2, (2) 
jk = dk <e ''bt p^pf^ - p^^-^ K . (3) 



The monomer density pi can be obtained from the conservation identity ([[]) that relates it 
to the other cluster densities. In these equations, pk = dpk/dt and Dj-Uk = ±[uk±i — Uk] are 
finite differences, t, dk and jk are nondimensional. t and dk are related to the dimensional 
time t* and decay coefficient dl as follows 

t = Qf, 4 = f. (4) 

Here the factor Q has units of frequency, it depends on the particular model we choose for 
dk, and will be determined later, jk is the net rate of creation of a + 1 cluster from a 
k cluster (the flux in size space), given by the mass action law. In we have made the 
detailed balance assumption to relate the kinetic coefficient for monomer aggregation to that 
of decay of a (/c + 1) cluster, dk- This implies that the equilibrium size distribution solving 
jk = has the form 



In and ©, is the binding energy of a cluster, required to separate it into its 
monomer components. Then the total energy measured with respect to a configuration in 



which all clusters are monomers is —J2k=2Pk^k- For spherical aggregates, 

Sk = ({k - l)a - ^a(A;§ - 1)) kBT. (6) 

This formula holds for /c 3> 1, but we shall use it for all k > 1. aksT is the monomer- 
monomer bonding energy which, in the case of precipitation of crystals from a solution 
or segregation by coarsening of binary alloys, r nay depend on the particle density p (vol- 
ume fraction) through some empirical formulas In Eq. a = 27s(47rt>^/3)^/(A;^T), 
where 7^ and v = V/M are the interfacial free energy per unit area (surface tension) and 
the molecular volume, respectively. Note that a and a are both dimensionless. The cor- 
rection 3o'kBT/2 in (jHI) ensures that ei = 0, and it improves the agreement between the 
nucleation rate obtained from the BDE and experiments Q|. More precise atomic models 
were proposed by Penrose et al [l^ . 

Eqs. (H)) , (0) , , © and a given discrete diffusivity dk form a closed system of equations 
that we can solve for an appropriate initial condition. If initially only monomers are present, 
we have pi(0) = p, and Pfc(O) = for k > 2. Before we obtain formulas for the kinetic 
coefficient dk, we shall recall the more salient features of the equilibrium size distribution. 

A. Equilibrium size distribution 

The equilibrium distribution (0) satisfies jk = and it can be written as 

Pfc = Ple-^^ (7) 
gk = -{k - 1) Inpi - ^ = ^a{kl - 1) - {k - 1) In (e» , (8) 

where Qk is the activation energy, equivalently given by 

3 2 

gk = (^k - {k - 1) (f, cTfc = -o-(A;3 - 1), (fc > 1), (9) 
¥P = ln(e». (10) 

Here ai = = gi. Assuming k ^ 1, gk achieves its global maximum g^ = ak'^^^/2 + 
ak^^^^ — 3cr/2 at the critical size 

* = (11) 
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FIG. 1: (a) Scaled activation energy gk/gm as a function of the scaled size k/kc- (b) Scaled 
dimensionless density r = />e" as a function of the scaled dimensionless monomer density ri = 
Pie" for the equilibrium distribution (solid line). Data correspond to liquid iron at maximum 
undercooling (dot-dashed line), whereas for disilicate glass, /O ~ pi (solid line). 



Eq. (jni) can be rewritten as 



gk/gm as a function of k/kc is depicted in Fig. Ufa). 

Rewriting the flux (jS)) in the BDEs in terms of the activation energy, we obtain 



3k 



4 |(e-^+9. _i) p^-D+pk). (13) 



Eq. (j2I) is a spatially discrete Smoluchowski equation with diffusion coefficient and drift 
velocity 



Vk 



dk (e-^+^'= - l) . (14) 



Notice that < for an activation energy that increases with k and f ^ > for decreasing 
Qk. Hence, gk indicates how the discrete advection Vk transports the clusters in size space: 
subcritical clusters shrink as time elapses while supercritical clusters grow with time. 
For the equilibrium densities (|7j), the conservation identity (P) becomes 



e^p=f2k {e'^pif e""^ = ^ A; e^'^-'^K (15) 



k=l k=l 



This series converges for e"pi = e*^ < 1 < 0), and diverges for e"pi > 1 > 0). At the 
critical micelle concentration (CMC), pi = e~° {ip = 0), we obtain the critical density above 
which equilibrium is no longer possible, 

oo 

e''p^ = l + Y,ke-''K (16) 

k=2 

For p > pc, the BD kinetic equations predict phase segregation, i.e., indefinite growth of 
ever larger clusters. 

B. The controlling parameters 

The simplest nucleation problem consists of solving the BD equations ((H), (j2I) and (fT!^ . 
with dimensionless activation energy gk = cTk — {k — l)(p, discrete diffusivity dk (to be chosen 
later) and initial conditions 

pi(0)=p, p2(0) = p3(0) = ... = 0. (17) 

The only parameters left in this initial value problem are p and a. p controls the long-time 
behavior of the BDE: If p < Pc given by ()16|) . pk{t) approach their equilibrium values Q, 
with monomer density pi that solves Equation (|T|): 

oo 

pe" = /(pie"; a)^J2 Kp^e'^f e-'^K (18) 

k=l 

The graph of this function is either the solid line or the dashed line in Fig. ^b). If p > pc, 
cluster sizes grow indefinitely whereas their density becomes small. Thus there remains a 
residual monomer concentration whose density pie° — 1 as t — oo. Summarizing, the union 
of solid or dashed lines in Fig. ^b) and the vertical line pie" = 1 for p > Pc represents the 
long-time limit of the monomer concentration as a function of p. 

Let us identify the controlling parameters p and cr in a physical system undergoing homo- 
geneous nucleation. A good experimental example is the transformation of certain silicate 
glasses to crystals (devitrification) In particular, abundant data exist for lithium dis- 
ilicate and we have compiled in Table I appropriate values of parameters characterizing 
nucleation j^. In disilicate, the free energy per molecule of the crystal phase in the activa- 
tion energy (jH)) is proportional to the undercooling 

_ Ag/(r„-r) 
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FIG. 2: (a) Comparison of s„(i) evaluated (at different times) from the numerical solution of the 
discrete equations to the asymptotic result ((51]) (solid line), (b) K{T) calculated from Eq. 
H41|) with K(0) = (solid line) is compared to the numerically obtained position of the wave front. 
Data correspond to disilicate glass at 820 K. All variables are written in dimensionless units. 



where is the melting temperature, ASj is the molar entropy of fusion and Na is Avo- 
gadro's number. The dimensionless density p = e'^^^^ can be extracted from Eq. (jl9p as ex- 
plained in Section lTTTl In energy units, the activation free energy is ksTgi^ = 7547ra^— fc^Tyj/c, 
where a is the radius of a spherical k cluster. From the expression for the volume of this 
cluster, kv = 4'7ra^/3 (v is the molecular volume), we obtain a = [3v/ {A7i)Y^^k^^^ , and 
therefore 

kBT (g,-^ + ^]= 7,(47r)i/3(3t;)^/3fc2/3 _ ^g^^j.^ _ T)k/NA. (20) 



2 

Comparing (pH) with yields a = {327iv^/3Y^'^js/{kBT), and the critical size 

~[ 3 ) ASjiT^-TY ^''^ 

The other parameters in Table I will be used later to model the discrete diffusivity in the 
BDE. We observe that the critical size increases with temperature: kc = 18 at 703 K and 
kc = 34 at 820 K. For other materials, such as undercooled liquid metals, critical sizes can 
be rather large: liquid iron at maximum und^ercooling has kc = 494, whereas kc = 2253 for 
liquid rutenium at maximum undercooling 



Parameter 


Symbol 


Value 


Melting temperature 


T 


1300 K 


Entropy of fusion 


ASf 


40 J mol-i K-i 


Surface tension 


7s 


0.15 J/m2 


Preexponential diffusivity 


Do 


2 xlO^ m^ s-^ 


Activation energy for diffusion 


Q 


440 kJ/mol 


Molecular volume 


V 


lQ-28 ^3 


TF time scale (703K) 


TF 


0.613 hours 


Heat capacity per unit volume 


PmC 


10^ J m-3 K-i 


Thermal conductivity (703K) 




3.96 X 10-^8 J m-i s-^ K'^ 


Thermally-driven-growth time scale (703K) 


TDG 


6.196 hours 


Critical size (703K) 


kc 


18 


Undercooling (703K) 


<P 


4.087 


Dimensionless surface tension (703K) 


a 


10.74 


Dimensionless free energy barrier (703K) 


n — £i-2/3 _ Sct _|_ 


25.177 


TF time scale (820K) 


TF 


0.0478 s 


Thermal conductivity (820K) 




1.84 X 10-^^ J m-^ s"^ K"^ 


Thermally-driven-growth time scale (820K) 


TDG 


0.48 s 


Critical size (820K) 


kc 


34 


Undercoohng (820K) 




2.817 


Dimensionless surface tension (820K) 


a 


9.207 


Dim('usi()ul(^ss free energy l)arri(M' ((S20K) 




38.181 



TABLE I: Data for lithium disilicate glass 



C. Equivalent Becker-Doring system 

As they stand, the BDE are rather stiff and hard to solve numerically. For example, at 
equilibrium. Table I indicates that Pk^/pi = e'^"""^ fa g-25.2-4.1 ~ 2 x 10"^^ for disilicate 
glass at 703 K, and pud Pi ~ 1-6 x 10""^^ for disihcate glass at 820 K. This motivates the 
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following change of variable 

Pfc = Pie-^'=Sfc = e-V'^-'^'=Sfc, (22) 

according to (fTUj) . Note that = 1 in equilibrium. Since gi = 0, this equation implies 

si = 1, (23) 

for all t. For the initial condition ()17|). e'^*^'^-'"" = pi(0) = p, and the conservation identity 
((T)) becomes 

oo 

k=2 

in which we have used fl22p . In terms of the Sk, the flux can be written as 

jk = dkexp[{k + l)ip - ak+i] {sk - Sk+i), (25) 
and the BDE (0) and ^ become 

s'k + Uk{sk+i - Sk) = -kipSk + dk-i {sk-i - 2sk + Sk+i), (26) 

for k > 2. Here, 

Uk = dk-i - dke'^-^'+'^K (27) 

The term Uk D+Sk in Eq. represents discrete advection, with a drift velocity Uk = 
—Vk + {dk-i — dk) ~ — "Wa;, which is essentially minus the drift velocity in the original BDE 
for k ^ 1. Thus, the advection in Eq. (j26p climbs up the activation energy barrier, from 
small values of Qk to large ones. 

In summary, the transformed nucleation initial-boundary value problem consists of the 
balance equations (j^Bj) . the particle conservation equation ()24|) . the boundary condition ^I^ . 
Si = 1, and initial conditions 5^(0) = for all k >2. Its solution gives ip{t) and Sk{t) for all 
> 2 and all t > 0. 

D. Stationary solution 

The stationary solution of the BDE has a flux independent of cluster size, so that jk = 
4exp[(/c + l)v9-(TA;+i] {sk-Sk+i) = j, from which {sk+i-Sk) = -j exp[ak+i- {k + l)(p]/dk, 
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and therefore 

,^^,_;ji'EEb±i^M2M, (28) 

1=1 

for k >2. Since Sqo = 0, j can be obtained from this expression in terms of an infinite series 

1 

^ ~ E£iexp[(Ti+i-(/ + l)v9-lnrfi]- ^^^^ 



Substituting back this expression into (|2H|) . we obtain 



^ . _ Eti^ exp[a;+i - (/ + - Indi] . 
Ez=iexp[a,+i-(/ + l)^-lnd,]- ^ ^ 



Then, p^ = p^e ^''Sk- 



E. TurnbuU-Fisher discrete difFusivity 

To solve the BDE, we need to estabhsh reasonable models of the kinetic coefficient dk 
(discrete diffusivity) for the decay of the {k + 1) cluster. A classical formula due to TurnbuU 
and Fisher 3] applies to spherical clusters whose growth is limited by the reaction rate at 
their boundary: dl is the product of the number of active sites on the aggregate times the 
molecular jump rate 



s the pr 

m 



an „2/3 „,2/3 Q/{RT) 

dl = A k^/' e^+^*/2^ = n e^+^^Z^ Q-^ = !— = . (31) 

Here D = Dqc^'^^^^'^^ is the diffusion coefficient in the liquid, Q is the activation energy 
for diffusion (see Table T), R = ksN^ is the gas constant and A = f^^^ {v is the molecular 
volume). If we nondimensionalize time as in Eq. (0)) with this definition of ^2, we obtain 

4 = fc2/3gDW2_ (32) 



III. ASYMPTOTIC THEORY OF TRANSIENT HOMOGENEOUS NUCLE- 
ATION WITH THE TURNBULL-FISHER DIFFUSIVITY 

In this section, we shall interpret the numerical solutions shown in Figures |21 to |3] by 
using singular perturbation methods. Our theory will be described using the TF discrete 
diffusivity ()31|) and compared to numerical solution of the BDE for the crystallization of 
disilicate glass at different undercoolings. 
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A. Initial transient 



Initially, pi(0) = p and there are no multiparticle aggregates. There is an initial transient 
stage during which dimers, trimers, etc. form at the expense of the monomers. This initial 
stage is characterized by the decay of the chemical driving force = « + In pi to a quasi- 
stationary value (p, given by Eq. ()19p in the case of disilicate glass, and the emergence of a 
continuum size distribution. Knowing this, we choose the initial chemical driving force v?(0) 
so that the quasistationary value (p given by Eq. Aiy^) is attained at the end of the initial 
stage. 

In materials such as disilicate glass at the temperatures we consider, the critical size is 
relatively small. Then v'(O) ~ ip, and the initial stage is very short. As the critical size 
increases (as in the case of undercooled liquid metals), ^^(O) may differ appreciably from (p, 
and the initial stage lasts longer. However, even in such cases, the duration of the initial 
stage, too, is negligible if we are interested in the overall duration of the transient stage 
to quasi- stationary nucleation. We shall show later that the duration of the initial stage 
compared to the duration of the overall transient is of order a very small quantity for 

materials with large critical sizes. 

B. Wave front advancing towards the cluster of critical size 

After the first era, clusters of increasing size are formed. For sufficiently small clusters, 
the continuum size distribution approaches the equilibrium distribution with ip = ip. This 
situation can be observed as an advancing wave front in the variable Sfc(t), satisfying ~ 1 
(equilibrium) behind the front and Sk ~ ahead of the front. This second era is described 
by Equations to (fTTj) with (p = ip and = 0. The critical sizes. 



for disilicate glass are relatively small, between 10 and 50, but they are large for undercooled 
liquid metals, generally between 100 and 1000. Hence we shall use as a small gauge parameter 




(33) 




(34) 
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Our asymptotic analysis will be carried out in the limit e ^ 0, and therefore kc = e — cxd. 
Then dk, Uk and au in Eqs. flH^ . and are smooth functions of /c > 0: 

u{k) = d{k - 1) - d{k) exp[(^ - a(A; + 1) + a{k)]. (36) 



^. Position of the wave front 

In the numerical solutions shown in Fig. Efa), the graphs of Sk vs. k at fixed time have 
clear inflection points at some k, where Sk ~ 1/2. The inflection point is taken as the 
position of the wave front. In the continuum model, the front position k = kf{t) is a smooth 
function which obeys 

kf = u{kf). (37) 

Scaling fc/ as 

fc; = - (38) 
(same scaling as kc = e~'^), the right hand side of (jSZI) becomes 

u{kf) = luiK)+Oie), (39) 

U{K) = 27^2/3 sinh (|(^"'/' - 1)) • (40) 
Eq. ()37|) can be rewritten as 

^ = t/(ir) ^ 2ir^/^ sinh (|(ir-V3 _ 1)) , (41) 

provided we deflne the slowly varying time scale T = et, and take the limit as e 0. 
Fig. Hfb) compares the position of the wave front calculated by solving (PT|) with K{0) = 
to the value obtained from the numerical solution of (j2Sl)- Note that the solution of ()41|) 
presents a time shift with respect to the numerical solution of the discrete model. This time 
shift reflects the breakdown of the continuum limit as K —>■ 0, due to discreteness, and also 
the transient in if{t) before it settles to (p. If the solution of Eq. (pT|l is forced to agree with 
the numerical K[T) when the latter is, say, 0.1, the comparison fares much better. 
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2. Shape of the wave front 



The leading edge of the wave front is a layer centered at K{T) in which Sk decreases from 
1 to as A; increases through it. The continuum representation of Sk in this layer is 



Sk = S{X,T] e) 



(42) 



where is 5* a smooth function of its arguments and X is the scaled displacement from the 
wave front location at A; = K/e^, i.e., 



(43) 

The scaling exponent p, presumably with < p < 3, is to be determined. The description 
()42|1 - ()43|1 should hold as e ^ with X fixed, so that the layer thickness scales as e~^. 
Substituting ^ into © yields. 



dS 



p-2 



dK dS 
dT dX 



K X 



+ «K + ^ [5(X + e^T;e)-5(X,T;e)] 



K X\ 



(44) 

Carrying out the straightforward expansion in powers of e, adopts the following asymp- 
totic form 



dS 



,2p-2 



p-2 



+ e 



(45) 



as e with X, K fixed. Here U{K) is given by (jiUI) . To obtain (j^ . we have used (jHHj) 
and PH): 

(K + e^~PX" 



u 



t-^U{K + e^-PX) + 0(e) = t-^U{K) + ei-Pf/'(K) X + o(e^- 



The dominant balance of diffusion and convection in yields 2p — 2 = 1, oi p = 3/2 
Hence (HHjl yields 



and the limit of as e ^ is 



dS ^^,.^^.^dS ^,r^.d^S 
^ + U'iK)X- = DiK)^, 



D(K) = lim 



e' = K'/\osh{^iK-'/'-r 



(46) 

(47) 
(48) 
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Had we carried out the same analysis for the ZFE, we would have found D{K) ~ d{e~^K) e. 
This would have resulted in a wider wave front and a longer time to ignition than those 
described below. 



3. Flux and wave front width 



Besides determining the shape of the wave front near its location, Eq. ()47|) yields the 
behavior of the flux (creation rate of clusters larger than k) near k = kf. If we substitute 
dna), (Sa) and (gUD into (ESI): 



= -/c^/^exp 



we obtain 
Here 



G{K) G'{K)X G'{K) G"{K) 



p3/2 



dS 
dX' 



G{K) ^ q> [^-K'/' - K 



(49) 



(50) 



is a scaled version of the activation energy 0. 

Since jk is proportional to dS/dX, it is convenient to differentiate ()47p with respect to 
X in order to obtain an equation for J = —dS/dX, 



(51) 



Notice that J is locally conserved, and the following integral conservation identity holds: 

dS 



1 = -[S] 



oo 
— oo 



dX 



:dX 



JdX. 



Eq. (fFT|) has Gaussian solutions satisfying 

1 



J(X,T) 



2AA(T) 



exp 



AAiT) 



(52) 



(53) 



which yields 



^(X,T) = -erfc 



X 



2JA{T) 



(54) 
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for the wave front profile. Inserting Eq. ()53|) in Eq. ()51|). we find tlie following equation for 
A{T) > 0: 

— - 2U\K)A = D{K). (55) 

Since K{T) is an increasing function, we can express A as a function of K. Inserting (j4ip 
in dnni), we get 

dA 2U'{K) _ D{K) 

dK ~ imy m)' ^ ^ 

Direct integration of this equation yields 

r dK 3i^4/3(l + gsinh2[£(^q^]) 
A = qU' + U' 0^^ = ^ \, ^ ^, (57) 



G"{K) 1 
2 ^ AA 



. (58) 



f/3 1^ 

in which q is an arbitrary constant. 

After insertion of ()53|). the flux ()49|) becomes 

;^2/3 g3^/(2.) r G'{K)X G'iK) 

3^ n I T \ 3 W2 ^ 

Here K = K{T) and A = A{T) are found by solving the differential equations (j41|) and (|55p 
with initial conditions K{0) = and A{0) = 3e^/(2(^), respectively. We have to set g = 
in Eq. (j57|) for A would become exponentially large as K = ^ otherwise. As T ^ oo, 
K ^1 and A^ -D{1)/[2U'{1)]. The definitions ^ and (gHl) of U{K) and D{K) imply 
f/'(l) = -^/3, £)(!) = 1. Hence, A 3/(2^) as T ^ oo, or as ^ 1 in Eq. §7^. The 
definition (jSOI) ofG{K) implies G(l) = <^/2, G'(l) = 0, and G"(l) = -(^/3. Hence, the limit 
as T ^ oo of the creation rate (j58|) is 



Notice that the terms proportional to X and have disappeared from this expression and 
therefore j^o is asymptotically uniform for X = 0(1). Eq. dSHl) is the classical quasi-steady 



nucleation rate of supercritical clusters due to Zeldovich 
from the stationary flux in the limit as e — > 0. 



17| , and it can be directly obtained 



C. The nucleation rate of supercritical clusters 

Let us now study the transient creation rate, in which j = j^^ increases from to the 
steady Zeldovich value (|^. As we have just seen, our theory predicts that the wave front 
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FIG. 3: (a) Evolution of the dimensionless flux at critical size j{t), and (b) number of clusters 
surpassing critical size Nc{t) as a function of dimensionless time for disilicate glass at 820K, kc = 34. 
Solid lines correspond to numerical results, dashed lines to the approximation given by Eq. ()58() . 
dot-dashed lines to the linearization approximation 1)67(1 and dotted lines to the approximation 
((C8() corresponding to linearizing the equations for K{T) and A{T) as in Appendix O 



profile is given by where K(T) and A(T) are solutions of Eqs. dTTp and / f33)) . re- 

spectively. The flux of clusters with sizes larger than k is then given by Eq. I[58^} . Setting 
k = kc = (critical size) and X = (1 — K(T)) / e^^"^ in this equation, we obtain the 
nucleation rate predicted by our theory, j{t). Its integral over time yields the number of su- 
percritical clusters, Nc{t). We shall consider now a different and more explicit approximation 
of these results. 

1. Linearization of the wave front speed about the critical size 

Let us fix k = kc = (critical size) in the definition of X: 



We now set X = K, in ()58)1 and perform the limit as e — > with k, fixed 



(60) 

. The result is 

(61) 
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provided we use the limiting stationary value (4A) ^ = —G"{l)/2. 



The transient turns on when k = [1 — 
the wave front equation (^01) yields 

dK 

It 

as e ^ 0. The solution of this equation is 



K)/e 



3/2 



0(1). Since [/(I -£3/2^; 



~ e 



3/2 ~ 



r 



-1 



-(pe. 



(62) 

(63) 
(64) 



It is convenient to choose km as the value of k, at which the flux j reaches its inflection 
point. Then we may consider that the wave front has ignited the nucleation of supercritical 
clusters. Straightforward use of Eqs. and (jU^ shows that 



Km 



(65) 



Moreover, Tm = is the time to ignition, at which the wave front K{T) reaches the value 
K = 1 — e^^'^KM- From (jHj), we obtain 



^00 ~l~ 



+ 



2<pe 



In 



'^(1-e 



3\2^ 



6e" 



+ 



K -1 



dK 



(66) 



37^2/3 sinh [f(ir-V3 _ 

where too is the duration of the initial stage. We could have expanded the integral in this 
expression, but Eq. is better suited for numerical calculation. The nucleation rate is 
found by inserting Eq. (jUHj) in (pT|) : 



J ~ joo exp 



'{t-tM)/T 



(67) 



in which a term of order e^/2 j^as been ignored in the exponential. 

Integrating j{t) over time, we find the number of supercritical clusters as a function of 
time. In the limit as t — > 00, this number is Nc{t) ~ joo (^ ~ ^)) where the time lag 6 is 
approximately given by 6* = + ^7 + rEi^e^^"/'^), or 



— too + 



+ 



In 



'^(1 



^3\2^ 



6e3 



7 



rEi(e*"/^) 



37^2/3 sinh f(ir-V3_ 1^ 



+ 



K 



dK 



(68) 
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FIG. 4: Same as in Fig. |21for disilicate glass at 703K, kc = 18. 

where 7 = 0.577215... is Euler's constant and Ei{x) is an exponential integral, see the 
derivation in Appendix El The time lag 6 can be directly compared to experimental values 

a. 

2. Comparison between different approximations 

Fig. Of a) compares j{f) calculated from the numerical solution of the BDE for devitrifi- 
cation of disilicate glass at 820 K, from (jU7j) and ()66|) with = 0, and from Eq. ()58|) with 
X = (1 - K(T))/e3/2. We find that the more precise expression, Eq. (j3H|) . captures better 
the width and location of the transition region between j = and j = joo, as compared with 
the simple approximation given by Eqs. (jHTj) and Both approximations present a small 
overshoot and yield a smaller time lag 6 than that obtained from the numerical solution 
of the BDE. The overshoot decreases as the critical size increases. Another approximation 
consists of linearizing the equations for K{T) and A{T) about the critical size K = 1 as 



suggested in Ref. 2M further explained in Appendix O This latter approximation is 
the worst one. This is not surprising as such approximation provides the same result for 
both the discrete BDE and the continuum ZFE. 
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For disilicate glass at a lower temperature of 703 K, the critical size is smaller and our 
approximations deviate more from the numerical solution of the BDE, as shown in Figure 
m^a). Integrating j{T) over time, we find the number of supercritical clusters as a function 
of time, Nc{t), which is depicted in Figs. E^b) and^J^b). At 703 K, the numerical solution 
of the BDE with the TF diffusivity yields a time lag 6 = 2.6. This value is close to those 
provided by the linearization approximation, 6 = 2.2, and by Eq. 6 = 2.3. Thus these 
analytical approximations to the numerical solution are reasonably good even for a relatively 
small critical size. However, 6 = 2.6 gives 1.6 hours according to Table I, whereas the 
experimentally measured time lag is about 50 hours, cf. Fig. 5 of Ref. ^. This discrepancy 
is due to having used the TF discrete diffusivity, which yields an excessively small time unit, 
as shown in Table I. 

IV. TEMPERATURE DRIVEN GROWTH OF THE NUCLEUS 

To improve agreement with experiments, we need a discrete diffusivity different from the 
TF one. We shall no longer assume that cluster size changes due to the activated transfer of 
a monomer through the cluster surface as in the TF theory. Instead, we shall assume that 
the discrete diffusivity in the BDE agrees with an adiabatic temperature driven growth of the 
nucleus. This yields a different formula for dk which, presumably, is not physically justified 
for very small cluster sizes. Nevertheless, the numerical solution of the BDE corresponding 
to thermally driven growth provides a time lag which is much closer to the experimentally 
measured value for disilicate glass than the TF diffusivity. 

A. Discrete diffusivity 

Let us assume that there is a nonuniform temperature field about a spherical crystal of 
radius a {k cluster) that is growing at the expense of the surrounding glass. Eq. (jTT)) shows 
that a nucleus of critical size grows if (T^ — T) decreases. The same equation yields (T^ — T) 
at the surface of a critical nucleus with k monomers: 



At the surface of the crystal, T = T^,{k), whereas far from it there is a smaller temperature, 
T = Too- Heat transfer from crystal to glass, —pmCK,4:na'^dT{a,t)/dr {k, c and pm are 




(69) 
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thermal diffusivity, specific heat and mass density, respectively), should equal the increase 
of energy due to crystal growth, TooASfNj^dk/dt*. We find 

dk pmCnNAdT ^ 

J r=a 

The temperature T(r) is the solution of Laplace's equation with boundary conditions T = T^ 
at r = a, and T = Too infinitely far from the nucleus. The corresponding solution is 
T = Too + (7"* — Too)a/r, which, together with (ffn|l . yield 

Using Eq. and the expressions for the radius a and (pij) for the critical size /Cc, we obtain 
rfA: 2(67r2t;)V3p^cft:iV^(T„ - Too. 



As ^ oo, the flux (jlHj) in the BDE becomes ~ d*^ (e"^^*^/^^ — 1) (written in dimensional 
units), and therefore 



{e^-^^-'" -l) = dl e-^'-'" (e^ - 



Using Eq. (fTTj). this equation can be written as 



dt 

Comparing Eqs. (|7^ and (|7Hjl . we obtain 

1/3 



'^\r^dl (e^[MW^)^/^]_i). (73) 



dl = nk^/^- 



_ 2(67r^i;)V3p^c/t fcgiVj 
As before, we shall absorb the constant f2 in the definition of time according to Q, which 



yields the following value of the dimensionless discrete diffusivity: 

1/3 



4 = k^'''- 



1 _ eV'[i-(fccA)V3]- ('^^) 

Here we shall assume that the thermal diffusivity follows the same Arrhenius law as the 
diffusion coefficient in the liquid k = D^e^^^^^'^^'^^] see Kelton et al 0^. With this choice 
of discrete diffusivity, a numerical solution of the BDE yields a time lag of 6' = 46.5 hours 
at 703 K compared to 9 = 1.6 hours previously obtained using the TP diffusivity. The 
experimentally measured time lag is 50 hours, as shown in Fig. 5, page 94 of Kelton's review 
j^. Thus we feel justified in using our formula f7^ to solve the BDE for disilicate glass. 

22 



B. Asymptotic theory 



We have to repeat the arguments given in Section IIIII using the discrete diffusivity (fTI)) 
instead of the TF expression. One important difference is that time needs to be rescaled 
as T = eH instead of T = et. Here we shall also use the symbol T for the slow time scale, 
but remembering that Ttdg = ^'^t (for thermally driven growth) instead of Ttf = (TF 
diffusivity). When necessary, we shall add the labels TF or TDG to the corresponding 
variables. After the initial discrete stage, our asymptotic theory yields the following results 
for temperature driven cluster growth, applicable to devitrification of disilicate glass: 

• The wave front profile S{X,T), with X = [k — e^^K(T)]e^^'^ and T = e^t, is given by 



(76) 
(77) 



S{X,T) = (l/2)erfc(X/[2./A(T)]). The front location and its width solve 



^ - 2U'{K) A = D{K) ^ I (1 - K"^) coth 



-1/3 



with initial conditions -ft'(O) = e^, A{Q) = e^/2. The latter condition corresponds to 
g = in Eq. ()57|) . Then A ~ K/2 a.s K —>■ 0+, which yields the initial condition for 
A if K{0) = e^. In Eqs. ^ and (HZj), U{K) and D{K) are defined by 



UiK) 



lim 



e u 



K 



D{K) 



lim 



d 



K 



1 fK 
—u 



Instead of Eq. (|49|) . we get the following approximation for the flux near the wave 
front: 



3k 



ei/2(^(l-iri/3)e3^/(2^) 



[1 



exp 



G{K) G'{K)X 



-G\K) 



G"{K) 

2 ^ AA 



g3/2 
X^ 



(78) 



in which G{K) = Lp (^K"^^^ /2 — K) and Eq. [\b'6\ has been used. Inserting X 
[1 — K{T)]/ e^/"^ in this equation, we obtain the nucleation rate: 

3{T) U{K) G{K) + G\K){l + t^ - K) 



I 3 U{K) 
2A(p 1 - e-G'(^) 
G"[K) 1 ■ 
2 ^ JA 



exp 



\2e3 
Kf\ 



'eo9 /3(y9 
— exp — 
Gtt ^ V 2e 



2e3 



(79) 
(80) 
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FIG. 5: (a) Evolution of the dimensionless flux at critical size j{t), and (b) number of clusters 
surpassing critical size Nc{t) as a function of time (in dimensionless units) for disilicate glass at 
820K, kc = 34. Solid lines correspond to numerical results, dashed lines to the approximation given 
by Eqs. CZni) to CZHl), and dot-dashed lines to Eas. (IHTl) - (IHHll. 



The simplest approximation for the nucleation rate and the time lag yields 

J = jooexp{-e-(*-*^^)/-}, (81) 

1 l-l-e^/^KM dK 



'TDG 



(^e2 y.3 1 _ i^l/3 

^ = er^F- (83) 

e = tM + b + Ei{e'^''^)]T 

To obtain these expressions, we have followed the same procedure as in the case of the 
TF diffusivity. In particular, Eqs. dHD) - ^ hold with T = eH, tn given by Eq. ^ 
and r given by Eq. 
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FIG. 6: Same as Figure [5] for disilicate glass at 703K, which has a critical size kc = 18. 

Figure El compares the numerical solution of the BDE (solid line) for devitrification of 
disilicate glass at 820 K (critical size kc = 34) with the more accurate asymptotic formulas: 
Eqs. dHU) - (j83|l (dot-dashed line) and its linearization about the critical size, Eqs. (|T6|l 
to (ffn|l (dashed line). Similarly, Fig. IHl corresponds to 703K. We observe that our two 
approximations, Eqs. (f7H|) - (|7n|) . and Eq. (jHTjl . describe quite accurately the numerical 
solution. Notice that our asymptotic formulas for thermally driven growth yield worse 
approximations to the numerical solution of the BDE than in the case of the TF diffusivity. 
The stationary nucleation rate is approximated less well by joo in the case of thermally 
driven growth because of the avoidable singularity of dk at the integer kc, which is slightly 
different from e~^. 

V. DISCUSION 

In this paper, we have studied the case of phase segregation resulting when p > pc- 
Previously, other authors had carried out asymptotic studies of the BDE in the simpler case 
of subcritical density, p < pc, in which initial conditions of only monomers, or more general 
ones, evolve towards the equilibrium distribution. In many cases of polynomial growth for 
dk, equilibrium is reached via a wave front profile for Sk, which is similar to Eq. with 
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A oc K^, and K oc T^, for appropriate positive S and /i; see Ref. j25[ and references cited 
therein. This advancing and widening wave front leaves in its wake the equihbrium size 
distribution. 

In the more complex case of phase segregation and indefinite aggregate growth considered 
in this paper, a quasicontinuum wave front of emerges after a short transient which is 
governed by the discrete BDE. After this, the leading edge of the wave front advances 
towards the critical size, and it slows down and stops there, leaving behind it a quasi- 
equilibrium state. The arrival of the wave front to the critical size marks the ignition 
of nucleation of supercritical clusters, which ends when the stationary Zeldovich rate is 
reached. Previous asymptotic theories have been derived for the continuum ZFE, not the 
discrete BDE, and thus their results systematically misrepresent two things: (i) the time 
lags for transient nucleation, as explained by Wu [19], and (ii) the width of the wave front 
and the time to ignition in the nucleation rate. The latter discrepancies occur because 
the diffusion coefficient appearing in the continuum equation for the wave front satisfies 
Dbde{K) = Dzfe{K) — U{K)/2, and therefore the time to ignition in the nucleation rate 
for the BDE is smaller than the corresponding one for the ZFE. 

Let us briefiy mention several existing asymptotic theories for the ZFE. Shneidman 



and Shi et al 



22| | Laplace transformed the continuum ZFE and matched a first stage of 



pure advection of clusters to a local expansion about the wave front when it is near its final 
position at the critical size. They obtained our simplest formula for the nucleation rate, 
Eq. (jnZj) with the same relaxation time, ttf or ttdg, except that their values for tM were 
different from This can be expected from Wu's arguments about approximating the 

discrete BDE by the continuum ZFE Q]; see the systematic shift of approximations of the 
ZFE with respect to numerical solutions of the BDE in Fig. 20 of Ref. . Trinkaus and 
Yoo 2o| studied a ZFE with a drift term linearized about the critical size (parabolic barrier) 
as an approximation to the full ZFE. Their results are comparable to those found bv means 
of the Laplace transform and matched asymptotic expansions; see Wu's review AH 
these authors obtained a transition region for the nucleation rate j{t) that was wider than 
observed in the numerical solution of the BDE. Several authors also found a nucleation rate 



for supercritical clusters that did not tend to joo as t ^ oo if A; 7^ j20|, |22, UM ; which is 



often called the asymptotics catastrophe j26|]. Our theory is free from this deficiency: Eq. 
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()B5|) in Appendix El provides the flux at k > kc using the TF diffusivity 



Joe exp 



-XoV2<^/3e 



-(i-iM)/(2r) 



e-^"'*"*«''^ (85) 



in which Xq = e"!"^ (k — kc). Notice that jk ~ joo as t oo, even after making our 
simplest approximation: hnearization of the wave front about the critical size. To get 
rid of the asymptotics catastrophe, Maksimov et al |3| assumed that S{X,T) = (1/2) 
erfc{[v4e-*/(2^) + B{X)]/ Jl^~Ce^}, in which the new function B{X) obeyed an ad hoc 



self-consistent equation that ensured jk ~ joo as t ^ oo even if 7^ kc- Note that if we 



use (j63|) for X = k, and the linearization approximation for A as in Appendix 
the previous formula for S with ( = I, A = e^^'^^'^'^^ and 5 = 0. Shneidman 



I we obtain 
2711 criticized 



Maksimov et al's result and extended his earlier asymptotic formula for the nucleation 

n 

rate ,23] to non-critical sizes. The previous criticism of using approximations to the ZFE 
instead of approximations to the discrete BDE apply to these works. Our more precise 
approximation using Eq. (j58|l plus the exact equations for the wave front location and its 
instantaneous width improve upon these approximations and perform better for materials 
with large critical sizes. 

The time lag obtained from the numerical solution of the BDE with the TF diffusivity (or 
from our asymptotic approximations using it) is too small as compared with experimental 
results (about thirty times smaller for disilicate at 703 K). The TF discrete diffusivity yields 
an excessively small time unit, as shown in Table I. We have greatly improved the agreement 
of theory and experiments by using a different formula for the discrete diffusivity, which is 
found by imposing that the growth rate of a critical nucleus resulting from the BDE be 
the same as obtained by heat transfer. In this case, our asymptotic approximations have a 
slightly different scaling of time and different expressions for U{K) and D[K). 
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APPENDIX A: GENERAL SOLUTION OF EQ. dSTJ 



It is convenient to rewrite this equation in terms of the variables K and X, as 

dJ U'{K)d{XJ) _ D{K) d'J 

dK u{K) dx u{K) ' y ) 

to be solved with the homogeneous boundary condition 

[' + U{K) ) U{K) dX - " ^^^^ 

(at X = Xin = e^/^ _/fg-3/2^ corresponding to A; = 1 in the definition of X), and with initial 
condition J(X, Kq) = —dSQ{X)/dX. The boundary condition is obtained by differentiating 

^(X,„,T) = 1 (A3) 

with respect to T and then using the definition of J and Eqs. ()41|) and (|Trj). The solution 
of the initial-boundary value problem is 



J(X, K) = - r G(X, Xo, K,) ^^^A dXo, (A4) 



dSoiXo) 



where the Green's function G{X, K; Xq, Kq) satisfies Eq. ()A1|) with initial condition 
G{X, Ko+; Xq, Kq) = 6{X — Xq) and the same homogeneous boundary condition as J at 
X = Xin- In a simple application of the method of images, the Green's function for this 
BVP can be written in terms of the Green's function Goo{X, K; Xq, Kq) for the infinite real 
line X, as 

G{X,K;Xq,Kq) = G^{X,K;Xq,Ko) 

+ c{Xq; K, Kq) G'oo(X, K- 2XM/U - Xq, i^o), (A5) 
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Vo J 




+ 1 


(^.-3/2+^' 


in 

U 


^2UB 




Uo ) 
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- 1 




U 


^2UB 



c{Xq;K,Kq) = ^ , (A6) 



Now, Goo{X, K; Xq, Kq) can be calculated by first writing an equation for the Fourier trans- 
form Goo(^, Xo, -K'o) = J^^e^^^Goo{X, K; Xq, KQ)dX . Such an equation is a first-order 
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quasilinear hyperbolic equation that can be solved by the method of characteristics. The 
result is that Goo is Gaussian in ^. Inverting the Fourier transform, we obtain 



G^iX,K;Xo,Ko) = , (A8) 

U^A7tB{K,Ko) 

Given the initial condition So{Xq) = H{2e^/'^ - Koe'^^"^ - Xq) (pure monomers), Eq. (jX4|l 
yields 

p 4172 B(K.Kq) 

J(X, K) = GiX, K; 2^1^ - e-^'^K^, Ko) , (A9) 



which is Eq. fl53|l . up to exponentially small terms. Here Kq = e^, Uq = U{Kq) 
2e2sinh[(^/(2e) - <f/2] ~ e2g^/(2.)-<^/2^ jj ^ u{K), and K = K{T). 



APPENDIX B: CALCULATING THE TIME LAG 



The time Tm can be estimated from Eq. ()41|) with initial condition -ft'(O) = (pure 
monomers) as 



U{K) 



dK 



3 U'{l){K~l) 



+ 



dK. 



Bl) 



U{K) U'(l){K-l) 

where too is the duration of the initial discrete stage in the original time scale. After straight 
forward calculations, we obtain (jHEl)- 
The number of supercritical clusters is 



Xc ~ 



j(t)dt = j^lt+ j exp(-e 



-{t-tM)h 



I JO 

-{t-tM)/r 



1 



1 — exp y—e 



dt = T dx 

Jo X 



dt\ r^joo{t-9), (B2) 



(B3) 



where Ei{x) is an exponential integral and 7 = 0.577215 ... is Euler's constant [28j. Notice 
that tEi{x) ~ re-^/a; ~ 9e[(^(l - e=^)]-2e-^/(6^') < e < 1, as a; = e^^'^^ ~ <^/(6e^) > 
Thus we can ignore the exponential integral in ()B3jl . which then yields Eq. fIfjSj) . 



1 Q- 
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Sometimes it is interesting to calculate the creation rate of clusters of size k > kc. If k is 
close to critical size, we can write 



X = e^lHk- 



3/2 



{k- 



+ K. 



(B4) 



Thus X = Xq + K, with Xq = e'l'^ (k - e~^). Inserting X = Xq + k and K =\- c'I'^k in 
(I58|) . we obtain the creation rate of clusters of size k = + X^e^^^'^: 



Joo exp 



XoJ2^/3e 



-(i-iM)/{2r) 



-(.t-tj^j)/T 



(B5) 



Notice that jk ~ joo as t — * oo. Thus our asymptotic result for the flux over any cluster size 
is free from the asymptotics catastrophe several authors found that their expressions 
for the flux tend to J 7^ joo as t ^ oo if A; 7^ kc- These catastrophes are due to inappropriate 
assumptions they made in their derivations. 



APPENDIX C: LINEARIZATION OF THE EQUATIONS FOR K{T) AND A{T) 
ABOUT THE CRITICAL SIZE 

A possible approximation of the wave front equations consists of linearizing the equations 
for K{T) and A{T) about the critical size K = 1: 

^^U\1){K-1), (CI) 
^-2U'{1)A^D{1), (C2) 

with K{0) = and A{0) ~ (initial condition of pure monomers). Both for the TF discrete 
diffusivity and for temperature driven growth of the nucleus, we have 

U'{1) = -|, D{1) = 1. (C3) 



The solutions of Equations ()C1|) and ()C2|) can be written in terms of the time scale t as 

ir(t) ^ 1 - (1 -e3)e-*/(2r), (C4) 
A{t) ^ ^ (1 - e-*/-) , (C5) 

for the TF diffusivity, and the same formulas with ttdg instead of ttf for thermally driven 
growth. Near the critical size, these equations would give an explicit expression of the wave 
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front profile (j54p . with X = e^/'^k — e '^/'^K{T). The nucleation rate of supercritical clusters 
is then obtained from with X = (1 - K)/e^/'^, together with dCH) and JHSl- For the 
TF diffusivity, j = jk, is 



1 - {l-e^)e 



-t/(2r) 



2/3 



e-t/r 



exp < — 



+ 



g3\g-t/(2T) 



-1/3 



(C6) 



2/3 



+ 



+ 



;i _ e3)2g-t/r 

6 



-t/r 



e3)e-t/(2r) 



-1/3 



-4/3 



(C7) 



Notice that the same results would have been obtained from the ZFE because the differ- 
ence between our D{K) and the corresponding one for the ZFE is U{K)/2, which vanishes 
at 7^ = 1. If we replace 1 instead of [1 — (1 



e^)e ^Z^^'^)] in the previous formulas, we find: 



J 



Jc 



^3\2 



-2t/T 



-t/r 



exp 



(C8) 



6e3 1 - e-*/^ j ■ 

Demo and Kozfsek's theory for the ZFE j2J| would yield Eq. (jC^Sj) for the nucleation rate 
once a couple of errors are corrected. They found t^k = 7ttf/5 instead of the correct 
relaxation time ttf, and an extra factor of 9 in the argument of the exponential in Eq. ()C8|) . 
Moreover, their exponential contains a factor (1 — e) instead of (1 — e'^). Demo and Kozfsek's 
figure 3 shows that their formulas do not improve as the cluster size increases, as one would 
expect of correct asymptotic expressions. Instead, they seem to optimize the nucleation rate 
of crystals in disilicate at kc = 27 (T=800K), as compared with numerical solutions. The 
earlier theory by Trinkaus and Yoo [2^ calculated the Green function for a time-dependent 
ZFE with a quadratic barrier and also used a linear equation for the position of the wave 
front. Thus their results are related to those in this Appendix. 
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